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Abstract 

We calculate the Isgur-Wise function by measuring the elastic scattering 
amplitude of a D meson in the quenched approximation on a 24^ x 48 lat- 
tice at /5 = 6.2, using an 0(a)-improved fermion action. Fitting the result- 
ing chirally-extrapolated Isgur-Wise function to Stech's relativistic-oscillator 
parametrization, we obtain a slope parameter = 1.2 . We then use 

-3 

this result, in conjunction with heavy-quark symmetry, to extract Vet from 
the experimentally measured B D*W differential decay width. We find 
\Vcb\\/TB f^-^Sps = 0.038 , where the first set of errors is due to experi- 
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mental uncertainties, while the second is due to the uncertainty in our lattice 
determination of p^. 
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Heavy-quark symmetry enables all the non-perturbative, strong-interaction physics for 
semi-leptonic B ^ D and D* decays to be parametrized in terms of a single universal 
function of lu = v ■ v', where v and v' are the four- velocities of the B and D mesons 
respectively [1,2]. This function, ^{lu), known as the Isgur-Wise function, is normalized at 
the zero-recoil point: ^(1) = 1 [2]. In order to determine the element Vcb of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix from experimental measurements of these semi-leptonic 
decays (which are made at > 1), it is necessary to know the Isgur-Wise function, and in 
particular its slope at the zero recoil point, ^'(1). In this letter we report on a lattice QCD 
calculation of the Isgur-Wise function, and on the corresponding determination of Vcb- An 
approach complementary to the one described here is currently being pursued by Mandula 
and Ogilvie [3] and by Aglietti [4]. These authors are attempting to formulate the heavy 
quark effective theory in Euclidean space, and to exploit this formulation for a numerical 
evaluation of the Isgur-Wise function. 

To obtain the Isgur-Wise function, we evaluate the elastic scattering matrix element 
{D(p')\cj'^c\D(p)) on the mass shell [5]. Because the electromagnetic current C7^c is con- 
served, this matrix element can be parametrized in terms of a single form factor: 

{D{p')\cYc\D{p)) = moiv + v'Y ]f\w) , (1) 

where p^'^ = mjjv^'^ and uj = v-v' is the four- velocity recoil. In the limit of exact heavy-quark 
symmetry this form factor is simply ^(lu). 

There are two sources of corrections to this simple result: 

h^\u)=[l + P^\u) + ^^\u)]au^). (2) 

The first correction, l3''^{uj), results from radiative corrections to the heavy-quark current. 
The second correction, j^''{uj), is due to higher-dimension operators with coefficients pro- 
portional to inverse powers of the charm quark mass. As defined in Eq. (2), ^{lu) is 
renormalization-group invariant and normalized to one at a; = 1 [6]. 

The radiative corrections can be evaluated analytically in a model-independent way: 
they are perturbative QCD corrections. To obtain these corrections we use Neubert's short 
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distance expansion of heavy-quark currents [6]. Neubert's result accounts for the full order 
as dependence of the heavy-quark current on the mass ratio, z, of the current's two heavy 
quarks. This is important because, in D — )• D transitions, z = 1 and order asz"' corrections, 
n = 1, . . . , oo, can be expected to be in the 10% range. 

Because Aqcd/ {'^^c) — 1/12, we would naively expect the corrections which are pro- 
portional to inverse powers of the charm quark mass to be of the order of 8%. These 
corrections are difficult to quantify because they involve the light degrees of freedom and 
are therefore non-perturbative. Luke's theorem [7], however, guarantees that there are no 
O (AQC£)/(2m,c)) corrections to h^\Lo) at zero recoil. Moreover, model estimates of these cor- 
rections appear to remain well below 3% over the range of experimentally accessible recoils 
[8]. We will neglect the C (AQC£)/(2mc)) corrections in extracting the Isgur-Wise function 
from h^^Lo) and will study the heavy-quark-mass dependence of h^'^Lo) in a later publication. 

Following Neubert [9] , we extract Vd from the experimentally measured differential decay 
rate for B — ^ D*li> decays. Using heavy-quark symmetry, we can express this differential 
decay rate in terms of ^{uj). In the hmit of zero lepton mass. 



where j5^'^ {uj) is the radiative correction corresponding to the form factor relevant for B — > 
D* transitions at zero recoil; P^'^iX) = —0.01 [6]. Moreover, up to corrections of order 
(AQc_D/(2mc,6))^, K(uj) = 1 at zero recoil (Luke's theorem [7]). Away from zero recoil, 
K{uj) contains l/rric^b and radiative corrections. Due to a fortunate cancellation, their sum 
remains small for all values of cu accessible in semi-leptonic decays [10]. Therefore we will 
neglect all non-perturbative effects in the coefficients which relate ^{lo) to the differential 
decay rate of Eq. (3), and set K{uj) = 1 for all lo. Then our lattice determination of ^{lo) 
enables us to extract Vet from the experimentally measured differential decay rate. 

We work in the quenched approximation on a 24^ x 48 lattice at (3 = 6.2, which cor- 
responds to an inverse lattice spacing = 2.73(5) GeV, as determined from the string 
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tension [11]. Our calculation is performed on sixty SU{3) gauge field configurations (for 
details see Ref. [11]). The mesons are composed of a propagating heavy quark with a mass 
approximately equal to that of the charm quark, and light antiquarks with masses around 
that of the strange quark. To reduce discretization errors, the quark propagators are calcu- 
lated using an C(a)-improved action [12]. This improvement is particularly important here 
since we are studying the propagation of quarks whose bare masses are around one third 
the inverse lattice spacing. Our statistical errors are calculated according to the bootstrap 
procedure described in Ref. [11]. 

The calculation of the matrix element (D(p') |c7'^c|D(p)) proceeds along lines which are 
now standard in the field of lattice computations of weak matrix elements. (For a recent 
review on this subject and references to the original literature see, for example, Ref. [13]). 
Thus, we calculate the ratio of three-point correlators, 

p. „^ ^ Ex,y e-^-^ -e-^P'-y(J^(t;,y)V-^(t,x) 4(0,0)) 
Ex,y(-/i.(^/,y)VO(t,x) 4(0,0)) 

where Jjj is a spatially-extended interpolating field for the D meson [14] and is the 
0(a)-improved version of the vector current cj'^c [15]. To evaluate these correlators, we 
use the standard source method [16]. We choose = 24 and symmetrize the correlators 
about that point using Euclidean time reversal [17]. We evaluate for three values of the 
Ught-quark mass (k; = 0.14144, 0.14226, 0.14262) which straddle the strange quark mass 
(given by Kg = 0.1419(1) [18]): one value of the heavy-quark mass approximating that of 
the charm quark (kc = 0.129 [19]); two values of the final D-meson momentum, p', and six 
values of the momentum, p = q + p', carried by the initial D-meson (these momenta are 
given in Table I). Data with momenta greater than (7r/12a)-\/2 are excluded because they 
have larger statistical and systematic uncertainties. To improve statistics we average over 
all equivalent momenta. 

Provided the three points in the correlators of Eq. (4) are sufficiently separated in time, 
the ground state contribution to the ratio dominates: 
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^ ' ^ ' ^^,i^-i^oo 2Eb-E^ Z|,(0) ^ ' ^ ^' ^^^^ ' 

(5) 

where we have used the fact that, for the continuum current, {D{0)\cj^c{0)\D{0)) = 2mD- 
En (E'jj) is the energy of the initial (final) D meson and the wavefunction factor, Znip^) = 
(0| j£)(0)|D(p)), is a function of the meson's momentum, because we use spatially-extended 
interpolating operators. 
We fit to 

R{t; p', q) = A\t- p', q) X e^Eo-E'^)t+iE'^-n.o)tf _ (g) 

This ratio becomes independent of t when both t and tf — t are sufficiently large, since the 
exponential factor in Eq. (6) explicitly cancels A^'s time dependence. We see a plateau in R 
about t = 12, and fit R(t; p', q) to a constant for t = 11, 12, 13. Multiplying this constant by 
suitable wavefunction and energy factors, we obtain {D{p')\cj^c\D(p)) . All wavefunction 
factors and energies are obtained from fits to two-point functions. Our results for ^{ui) are 
presented in Table I. 

The data for ki = 0.14144, the heaviest of our light quarks, which have the smallest 
statistical errors, are shown in Fig. 1. The solid curve is a two-parameter fit to s^p{uj), 
where ^p{lo) is Stech's relativistic-oscillator parametrization [8,10]: 

«»-^-p(-(V-l)^) (7) 

and = —(^^(1). The parameter s was added to absorb uncertainties in the overall nor- 
malization of our data through a common factor. We find = 1.5 and ,s = 0.95 with 
a x^/do/ = 1.0. Other parametrizations for ^(lu) give nearly identical results; for instance, 

H~ 2 -|- 1 

the pole ansatz of Ref. [10] yields = 1.4 and s = 0.95 with a x^/'iof = 1.0. The fact 
that s is not quite consistent with 1 indicates that there may be some small systematic un- 
certainty associated with our choice of normalization, or that the standard parametrizations 
for ^{uj) are not optimal. 
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Now, if we fit our data to ^p{lu) instead of s^p{uj) (an equally valid procedure, in principle, 
for determining p^), we find = 2.0'*'^ with a x^/do/ = 2.6 (dotted curve in Fig. 1). To 
accommodate the spread in values for given by our two procedures, we assign errors to p^ 
which encompass all values consistent with both procedures. These errors include systematic 
uncertainties, but only to the extent that the deviation of s from 1, in our first fit, is an 
indication of systematic errors. The central value we choose for p^ is the one given by our 
first fit since this fit is designed to absorb possible uncertainties in the overall normalization 
of our data. Thus, for ki = 0.14144 we quote p^ = 1.5 ^. 

In Fig. 2 we plot the results obtained from a covariant and linear extrapolation of our 
data for the three values of the light-quark mass to ki = Kcrit- These results correspond to 
a meson composed of a charm quark and a massless antiquark. Since ^{uj) is not thought 
to depend very strongly on the light-quark mass [20], we expect the slope parameter, p^, 
for this data to be close to that found when = 0.14144, i.e. when the light antiquark is 
slightly heavier than the strange quark. This is indeed Avhat we find. With a two parameter 
fit to s^p{uj), we get p^ = 1.2 and s = 0.94 with a jAoj = 0.9. Forcing s to be 1, 

we find p^ = 1.7^'^ with a x^/do/ = 1.6. Using the same procedure as the one used above 
to determine the errors and central value for p^ when ki = 0.14144, we obtain, as our best 
estimate for p^, when ki = Kcrit- 

P' = 1.2^_l. (8) 

Our result for p^ agrees with most other determinations of this parameter [9,21,22], apart 
from the sum-rule result of Ref. [23] which lies below our error bars. In particular, our result 
for p^ agrees with the very recent lattice result of Ref. [24] obtained with Wilson fermions, 
although the details and systematics of the two calculations are different. 

Having chosen a parametrization for the Isgur-Wise function (Eq. (7)), and having de- 
termined the slope parameter (Eq. (8)), we can now obtain Vcb- Setting K{(j) to one, as 
discussed after Eq. (3), we fit the decay rate of Eq. (3) to the experimental data. In Fig. 3 
we show a least-x^-fit to the new ARGUS data for B — > D*W decays [25]. The resulting 
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value of Vcb is 




1.48ps 



0.038 



+ 2 + 8 
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(9) 



with a x^/do/ = 1.1. The same fit to the weighted average of older CLEO and ARGUS 



set of errors is due to experimental uncertainties, while the second is due to the uncertainty 
in our lattice determination of p^. The B-meson lifetime used above is the central value of 
the lifetime, r^o = 1.48(10)ps, quoted in Ref. [27]. 

In this letter we have reported on a lattice computation of the Isgur-Wise function and 
on the corresponding determination of the CKM-matrix element, Vcb- It should be stressed 
that the lattice errors for in Eq. (8) and for Vcb in Eq. (9) include systematic errors only 
to the extent that the spread in values for given by the two procedures discussed after 
Eq. (7) is a measure of these errors. We have tried to minimize systematics by working with 
an improved action to reduce discretisation errors, and on a fairly large volume in the hope 
that finite size effects would be small. Nevertheless, it is important that our simulation be 
repeated on lattices with different lattice spacings and on lattices of different sizes in order 
to quantify more precisely these systematic effects. It should also be remembered that the 
computation was performed in the quenched approximation. 

Lattice simulations provide the opportunity for testing the validity of the Heavy Quark 
Effective Theory for charm quarks, by allowing one to check whether is independent 
of the mass of the heavy quark. They further enable one to check the extent to which 
the Isgur-Wise function obtained from pseudoscalar — )• vector transitions agrees with that 
obtained from pseudoscalar — ^ pseudoscalar transitions. Both these checks are currently 
under investigation. 



data [26] gives | ^6 1 ^Wl^^Sps = 0.036^^^^ with a xV^o/ = 0.6. In both 



cases, the first 
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FIGURES 

FIG. 1. The squares are our lattice measurements for the Isgur-Wise function, ^{lo), at the 
heaviest hght-quark mass (k; = 0.14144). The sohd curve is obtained by fitting our measurements 
to s^p(cj), where ^p{io) is Stech's relativistic-oscillator parametrization for the Isgur-Wise function 
(Eq. (7)). The dotted curve is obtained by fitting our data to 

FIG. 2. The octogons are our lattice results for the chirally-extrapolated Isgur-Wise function, 
i.e. ^{lo) for = Kcrit- The solid curve is obtained by fitting these results to s^p(tj), where Cp{io) 
is Stech's relativistic-oscillator parametrization for the Isgur-Wise function (Eq. (7)). The dotted 
curve is obtained by fitting our results to £,p{io)- 

FIG. 3. The best fit of |Fc6|^p(c(;) to experimental data, where is fixed to its value in Eq. 
(8). The experimental data (diamonds) are obtained from ARGUS's new results for B — )■ D*lu 
decays, assuming a B meson lifetime of 1.48 ps [27]. Also shown are our appropriately scaled, 
chirally-extrapolated results (octagons). 
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TABLES 

TABLE I. Results for the Isgur-Wise function, ^{lo). Only data for our heaviest light quark 
{ki = 0.14144) and for a massless light quark (k; = Kcrit= 0.14315(2) [18]) are displayed. The 
hopping parameter of the heavy quark is k/, = 0.129. p (p') is the momentum of the initial (final) 
pseudoscalar meson in lattice units. 
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